DIFFUSIVE PROPAGATION OF WAVE PACKETS IN A FLUCTUATING 

PERIODIC POTENTIAL 



EMAN HAMZA, YANG KANG, AND JEFFREY SCHENKER 

Abstract. We consider the evolution of a tight binding wave packet propagating in a 
fluctuating periodic potential. If the fluctuations stem from a stationary Markov process 
satisfying certain technical criteria, we show that the square amplitude of the wave packet 
after diffusive rescaling converges to a superposition of solutions of a heat equation. 



1. Introduction 

It is generally expected that wave packets evolving in a homogeneous random environment 
propagate diffusively over long time scales, unless recurrence effects are strong enough to 
induce Anderson localization. If furthermore the environment fluctuates in time, recurrence 
effects should be irrelevant, suggesting that diffusion is universal for wave motion in time 
dependent random systems. This idea was confirmed by Ovchinnikov and Erikman [5], 
who showed diffusion for a tight binding Schrodinger equation with white noise potentials. 
Pillet [6] considered a more general setting in which the potentials were Markov processes, 
but not necessarily white noise. He demonstrated the absence of binding and derived a 
Feynman-Kac formula. This Feynman-Kac formula was used by Tcheremchantsev [8] 
to show that position moments scale diffusively up to logarithmic corrections. Recently, 
two of us proved diffusion of wave packets and diffusive scaling [1] for a large family of 
Markov models, including those considered by Tcheremchantsev. For a recent discussion of 
the physics and physical applications of the tight binding Schrodinger equation with time 
dependent randomness we refer to 

The study of diffusion for disordered quantum systems, or "Quantum Brownian motion," 
has recently attracted the attention of a number of authors. Diffusion in the presence of a 
weak static random potential for a quantum particle on a lattice of dimension three or higher 
has been demonstrated only up to a finite time scale proportional to an inverse power of the 
disorder strength [U El E] • Frohlich, Pizzo and De Roeck have proved diffusion, for arbitrarily 
long times, for a quantum particle on a lattice weakly coupled to an array of independent 
heat baths [9J. In [9], it is mentioned that the method used therein also applies to a particle 
in a time dependent potential provided one has exponential decay of time correlations, such 
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as one has for the Markov potentials in [3]. However the proof in [9] relies on a polymer 
expansion which restricts the result to weak coupling (this is not the case in [3]). Another 
result closely related to our previous work [3] is a recent paper on diffusion starting from a 
quantum master equation in Lindblad form [TO] . 

This note and the aforementioned [61 [7J El [3] are concerned with the evolution of wave 
packets for the "tight binding Markov random Schrodinger equation:" 




where 

(1) L is a translation invariant hopping operator on £ 2 (Z d ), 

(2) v x : — > M are real valued functions on a probability space Q, 

(3) oj(t) is a Markov process on Q with an invariant probability measure //, and 

(4) v x (u) = Vo(a x (u))) where a x is a group of [i- measure preserving transformations of Q. 

(Formal definitions are given in section [2] below.) 

The potentials considered by Tcheremchantsev [7J E] were independent at different sites. 
However, this played no role in the analysis in (3]. Nonetheless, some non- degeneracy as- 
sumption is certainly needed as can be seen by considering the case v x = vq for all x, for 
which the effect of the random potential is only to multiply the wave function by a time 
dependent random phase. The technical condition employed in [3] was 

(1.2) inf IIS-^-OlL 2 ,_. >0, 



where B is the generator of the Markov process u(t). 

Our aim here is to consider a situation in which (11.21) is violated in a relatively strong 
way. Namely, we shall consider periodic potentials, v x+ ^ y = v x for all x, y with N some fixed 
number. Because the resulting system is periodic under translations by elements of NZ* d , 
there is a conserved "quasi-momentum." Our main result, in short, is that after taking into 
account conservation of quasi-momentum the motion of the wave packet is diffusive. More 
specifically, over long times one sees a superposition of diffusions: 

(1.3) lim y e "^ k "E(|^ Ti (i)| 2 ) = / e- iE ^=i D ^ )k ^m(p)dp, 

x£Z d l N 

where = [0,27r/iV) d , p \-> Aj(p) 

is a continuous function taking values in the positive 
definite matrices, independent of ipo, and 
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with A = [0, N) d n Z d . The quantity m(p) is the amplitude of the initial wave packet at 
quasi-momentum p — ipo denotes the Fourier transform of i[)q: 

(1.5) ft(k) = ^e-Vo(x), 

if ^ 6^n £ 2 . 

To understand the meaning of (11.31) . consider the following position space density 

(1.6) dRt(x) = J2 E (l^(e)| 2 ) 6{x - £)dx, 

a probability measure on IR d . (Here 5(x)dx is the Dirac measure with mass 1 at 0.) After 
taking inverse Fourier transforms of both sides, ( II. 3p shows 



;i.7) / <f>{x)dRrt(V¥x) > 

J Rd t^oo 



[ (f)(x) I - d — 1 - e >: •■ /; ■ !;,> "' '' ///(p)dp 

Jj d N (47rt)2^/det Aj(p) 



dx, 



for any test function on M. d which is, say, smooth and compactly supported. The function 
appearing as the integrand inside square brackets on the right hand side is the fundamental 
solution to an anisotripic diffusion equation, with diffusion matrix JDjj(p), 

(1-8) §«,(*) = E«^(4 

Thus (I1.3P can be understood as saying that the position space density dR t (x), after diffusive 
rescaling t rt and x t— > v^ 1 ^; converges in the weak* sense to 



1.9) dRrti^Fx) ^ 



u t (x; p)m(p)dp 



dx, 



where «t(x; p) satisfies (II. 8p with Uq(x; p)dx = <5(x)dx. That is over long time scales, after 
diffusive rescaling, the mean square amplitude breaks into components for each p, with each 
component propagating independently and according to a diffusion equation, which is to say 
a "super-position of diffusions." 

The result is stated formally in section |2] after we give the required assumptions. These 
assumptions are somewhat technical, so it may be useful to have a simple example in mind. 
Fix a function U : Z d — > R periodic under translations in A r Z d , that is, U(x — Ny) = U (x) 
for all x, y G Z d . Now let u(t) be a continuous time random walk on A = [0, N) d fl 7L d taken 
with periodic boundary conditions and with independent identically distributed exponential 
holding times at each step. The probability space is just A with the measure \i normalized 
counting measure. Take the potentials v x to be v x {uj) — U(x — uS) so that the Schrodinger 
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equation describes a particle in a "jiggling" periodic potential: 

(1.10) idMx) = Y, h (OMx-0 + u(x-u(t))Mx)- 

c 

Our result shows that (11.31) holds provided U has no smaller periods, i.e. that 

^2\U(x + y)-U(y)\^0, 16A and i/0. 

2/GA 

2. Statement of the main result: A superposition of diffusions 

2.1. Assumptions. Our main result is formulated with the following assumptions. (See [1] 
for a more detailed discussion of the framework.) 

Assumption 1. We are given a topological space f2, a Borel probability measure /i, and 
a Markov process on Q with right continuous paths for which /i is an invariant measure. 
Furthermore, we suppose that there is a representation of Z d , x i— >■ cr^, in terms of //-measure 
preserving maps <r x : f2 — > Q such that the paths of a x (u>(-)) have the same distribution as 
the paths of cu(-), for all x G Z d . 

We denote by E (•) expectation with respect to the paths of the Markov process with the 
initial condition cu(0) distributed according to /i. By the invariance of /i, we have 

(2.1) E(/(w(f)))= / f(a)M<*) 

Jn 

for any t and any / G L x (f2). Furthermore, the map given by 

(2.2) SJ(a) = E(f(u(0))\u(t) = a) 

defines a strongly continuous contraction semi-group on L 2 (Q). By the Lumer-Phillips theo- 
rem, St is generated by a maximally dissipative operator B with dense domain T>(B). Since 
Stl — 1 for all t, 51 = and is an eigenvalue of 5. Since B is dissipative, we also have 
that its numerical range lies in the right half plane. We suppose further that B is sectorial 
and satisfies a "spectral gap" condition: 

Assumption 2. There exist 7 < 00 and T > such that 

(2.3) \lm(f\Bf) L 2 (n) \ < 7 Re(/, Bf) L 2 {a) , 

and 

(2-4) Re</,£/) L2(n) >ivar(/) 

for all / G where Var(/) := /„ f 2 dfi - {f Q fdfi) 2 . 

The potential v x : — > R and hopping operator L are assumed to be translation invariant, 
and L should satisfy a non-degeneracy condition that precludes hopping only in a sub-lattice: 
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Assumption 3. The potential is given by Borel measurable bounded functions v x : Q — > R 
such that 

v x = v o a x . 

The hopping operator is given by 

Lip( x ) = ^h(x -y)ip(y), 
y 

where h(—x) = h(x)*, \x\ 2 \h(x)\ < oo, and for each non-zero vector k G IR d , there is 
some x G Z d such that h(x) 7^ and k • x 7^ 0. 

Finally, since we are concerned with periodic potentials, we suppose 

Assumption 4. There is N G N, N > 1, such that a^ x = Id for all x G Z d . Furthermore, 
we suppose that \\v x — v \\ Loa ^ > for all x G [0, N) d fl Z d , x 7^ 0. 

Remark. More generally, we might allow different periods in each of the coordinate directions: 
Ni, . . . N d such that a y = Id whenever y = (A^ax, . . . , N d ad) with a±, . . . , G Z. The result 
stated below holds also for this case with essentially the same proof. We choose to work 
with equal periods for notational clarity. 

Let A = [0, N) d n Z d , as above, and let x, y G A. Since v x — v y is mean zero, it is in the 
domain of B~ l . Furthermore, it follows from Assumption 4 that \\v x — Vy\\ L 2nm 7^ if x 7^ y, 
in which case — v y )\\ L2 ^ 7^ 0. Since A is finite, we conclude that there is x > 

such that 

(2-5) \\B~ l (v x -v y )\\ L2{u) >x, x,yeA,x^y. 

Eq. (12. 5p will play a key role in the proof below. 

2.2. Main result. Consider the density matrix 

(2.6) Pt(x,y) = ip t (x)ip t (y)* . 
It is well-known that pt(x,y) satisfies 

(2.7) d t p t (x,y) = -i^2h(()[p t (x - (,y) - p t (x,y + (v x (u(t)) -v y (u(t))) p t (x,y). 

C 

More generally, we may consider solutions to (12.71) with an initial condition 

(2.8) p e VM := {p:Z d x Z d — 7- C : p is the kernel of a non-negative definite, 

trace class operator on i 2 (Z d ) j . 

Recalling the notation = [0, 27r/N) d , we now state our theorem. 
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Theorem 1. The solution to ( 12 .TP initial condition po G D.M satisfies 
(2.9) lim Ve" i ^' x 'E(p T i(x,x)) = / e"' E ^ D ^ (p)k ^m(p)dp, 

where p i— >■ -Djj(p) zs a continuous function taking values in the positive- definite matrices 
and 

m(p) = ^/(*P), 
wt/i / i/ie Fourier transform of f(x) = J2 yeZ d Po(y + Nx,y). 

Remark. We have defined the function m in terms of the Fourier transform of /. Since / is 
not obviously summable or square summable, it is not immediately clear that m is indeed a 
function, rather than a distribution. However, in terms of the orthnormal eigenvectors ipj of 
Po and corresponding eigenvalues Xj, we have 

( 2 - 10 ) ™(p) = -^y E a, g |s (p + |c) I • 

. . 2 

Since Xj < oo and -0j G L 1 (T d ) we see that m(p) is an L 1 function of p. (The function 
f(x) can be expressed as 

(2.11) f(x) = tr Po S Nx 

where po is interpreted as a trace class operator and S^ x is the shift by Nx on £ 2 (Z d ), 
SNx^hj) = ^(y ~ Nx). It follows that / is positive definite: 

n 

(2.12) J]C*Ci/(^-^)>0 

for any finite collection of points x±,...,x n G Z d and any (Ci, • • • , Cn) £ C n . We conclude 
from Bochner's theorem that / is a non-negative measure of mass /(0) = tr p , and because 
lim^^oo /(x) = the measure has no point component. But, it is not immediately clear that 
/ is absolutely continuous with respect to Lebesuge measure so that m is a function. For 
this purpose (12.101) seems to be necessary.) 

3. Augmented space analysis 

In this section, we explain briefly the augmented space analysis, which was also employed 
in [H Section 3]. We begin with the following Feynman-Kac formula [6] 

(3.1) E(p t (x, y)) = (5 X ® 8 y <g) 1, e~ tL po ® 1>«, 

which relates M(p t (x, x)) to a matrix element of a contraction semigroup e~ tL on the aug- 
mented Hilbert space 

(3.2) U := L 2 {Z d x Z d x ft). 
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The operator L in (13. li is given by L := \K + iV + B, where 

(3.3) K*(x,y,oj) = J2 h (0[*(z-{,y,u)-V{x,y + {,co)}, 

c 

(3.4) ^*(2C,j/,w) = (^(w) -v„(a;))*(a;,j/,a;). 

The Markov generator B acts on H as a multiplication operator with respect to the first two 
coordinates: 

(3.5) B[p®f] = p® (Bf), p e t\Z d x z d ), f e L 2 (Sl). 

Our analysis, as in [I] , makes crucial use of the invariance of the generator L with respect 
to simultaneous translation of position and disorder. In the present context, we have a 
larger group of symmetries due to periodicity. Namely, the generator L and its constituents 
K, V, and B, commute with a group Q of unitary maps on H generated by the following 
transformations: 

(1) Simultaneous translation of position and disorder by an arbitrary element of Z d : 

S(V(x, y,u) = '9(x-£,y-£, cr^u), 

(2) Translation of the first position coordinate by an element of NZ d : 

Note that S^S^ = S$, so the group Q is isomorphic to Z d x Z d . We have chosen to use 
translation of the first position in the definition of however, since cr/v£ = Id, we have 

s n\ = SnzSH%£ G Q, where S^\^(x, y, u) = V(x t y - N£, u). 

Because of the invariance with respect to Q, L is partially diagonalized by the following 
generalized Fourier transform: 

(3.6) (x, u, k, p) = e^'-^-^Vix -t-N V , ap), 

a unitary map from L 2 (Z d x Z d x fi) L 2 (A x Q x Tf x T^). Thus we have, by flBTTD . 

N d f 

(3.7) J2e- ik - x E(p t (x,x)) = -— d dp(5o®l,e- iLk 'Pp ;k,p®l)L 2 (AxC), 
where 

(3.8) p ;k, P (x) = J2^< x - N ^- ik - y p (x-N V -y,-y), 

and Lk P := iK^ p + iV + B with 

(3.9) Vi>(x, u) = (v x (u) - vq(u))4)(x, lo), 
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and 

(3. 10) #k lP fe u) = ^2 HCy P< - C, - e" ik< ^(x - C, ^) 

C 

(In (13.101) we take "periodic boundary conditions," that is x — C on the right hand side is 
evaluated modulo N.) 

The transformed Feynmann-Kac formula (13.71) is the starting point for our proof of The- 
orem [Tj It reduces the study of the mean density in (12. 9p to the spectral analysis of the 
semi-group e~* Lk p for each fixed p and for k in a small neighborhood of 0. 

4. Spectral analysis of L kp and the proof of Theorem [TJ 

In this section inner products and norms are taken in the space L 2 (A x Q) unless otherwise 
indicated. We denote by Pq the orthogonal projection of L 2 (A x Q) onto the space "Ho = 
£ 2 (A) £g> {1} of "non-random" functions, 



(4.1) P V(x) = / *(x,u)dfi(u), 

JQ 

and by P^~ = (1 — Pq) the projection onto mean zero functions 

(4.2) H£ = U(x,u) 



J *(x,w)d/i(w) = o| . 



A preliminary observation is that 
(4.3) Z ,p<5 ® 1 = 

for all p. Thus, 5q <S> 1 is stationary under each semigroup e~ tL ° p . Eq. (14.31) can be seen 
easily from the explicit form for L p given above, but could also be derived from the fact 
that, for each y e Z d , 

J2^(pt(x + Ny,x)) 

X 

is constant in time. 

A key step toward proving Theorem [TJ is to observe that the remaining spectrum of L p 
is contained in a half plane with strictly positive real part. To see this, we make use of the 
block decomposition of L , p with respect to the direct sum 1~Lq © T^q ". 



(4.4) L 



p 



iP V 
iVP YK v + B + iP^VP£ 



(Note that K 0v and B both act trivially on T-L , while PqVPq = since J n (v x (u) — 
v (uj))d^(u) = 0.) 

We use (14. 4p to prove the following 
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Lemma 2. There is 5 > such that for all p G T%, 

(4.5) a(L 0>p ) = {0}U£+(p) 

where is a non- degenerate eigenvalue and £+(p) C {z : Re z > 6} . 

Proof. This is very close to jH Lemma 3]. The key new point is that we must see that 5 can 
be chosen independently of p. 

Because Re 5 > ^Pq, it follows from an argument using Schur complements that a point 
z with Re z < ^ is in a(L 0p ) if and only if z is in the spectrum of 

(4.6) T p (z) = P V(P ± L ^P ± - zy'VPo. 
However, given G ^ 2 (A), 

Re(0® l,r p (z)0® 1) 

, (P - L Zo,pi 3 1 - z) -1 ^ ® 1 , (Re B - Rez)(P ± Z ,pP ± - z)- x V> ® 1 



(4.7) 



(B- 1 P L ( J L Ol p - z) ^o ) 5 ® 1 



where the inverses are well defined because V"</> ® 1 G = ranPj 1 . Since ||P 1 -Po" L || — ^ 
it follows that 



(4i 



B- 1 P J -(L , p -z)P ± < 1 + T ||Ar ,p 



+ k 



However, ||i£k,p|| < 2 1| /i || oo for all k and p, so P 1 P^(L 0p — z)Pq is uniformly bounded and 
Re(0® l,r p (z)0® 1) 



(4.9) 

Finally, 
(4.10) 

where 
(4.11) 



I f~^z 



P"V0® l 



i + r(2||/i|| 00 + 2||v|| + ^|) 



P"V0® l 



E 



|p X K -^o)!!^^ > x 2 J2 



min IIP 1 (v x - v )\ 



L 2 (n) 



which is positive by Assumption 4. 
Thus, 



(4.12) 



l^-Rez 



l + T(2\\h\\ 00 + 2\\V\\ + \z\ 
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Since the right hand side is independent of p, the existence of a spectral gap 5 independent 
of p, as claimed, now follows from the sectoriality of B (Assumption 2, eq. ( 12. 3p ) as in the 
proof of [U Lemma 3], with the explicit estimate 

1 Y 2 

(4.13) 5>- - . □ 

T (2 + 7 + 4T||£|| 00 + 4T||V||) + \\V\\ V 

4.1. Analytic perturbation theory for L^p. We now hold p fixed and consider the 
spectrum of p for k close to 0. We write V for the gradient with respect to k and di for 
partial differentiation with respect to the i th coordinate of k. No derivatives with respect to 
p appear below. 

The key observation is that the spectral gap for L 0p is preserved in the spectrum of L k p 
for k sufficiently small. 

Lemma 3. Given e e {0,5), with 5 as in Lemma\^ there exists r such that if |k| < r then, 
for each p G T d N , 

(1) Lk, p has a single non- degenerate eigenvalue E p (k) with < Re_E p (k) < 5 — e, 

(2) The rest of the spectrum of L^ p is contained in the half plane {z : Re z > 5 — e}. 
Furthermore, -Ep(k) is C 2 in a neighborhood ofO, 

(4.14) E p (0) = 0, V£ p (0) = 0, 
and 

(4.15) didjE^O) = 2Re(d i K Oi p5 o 0l,[L Oi p]- 1 d J K OtP 5 o ®l) 

= 2Re x iyjW)Hy)(8[x] N ®h\rp(o)]~%y®i), 

where [x]n denotes the point in A equivalent to x modulo N and 

(4.16) T p (0) = P V(P L Lo, p P A -)- 1 VP . 
In particular, didjE p (0) is positive definite. 

Proof. These are essentially standard facts from analytic perturbation theory. The key point 
is that 



(4.17) 



7 k,p 



< c\k\. 



If the generators L^ p were self-adjoint or normal it would now follow that the spectrum 
moves by no more than a distance c|k| for k small. However, Lk )P need not be normal so we 
must argue more carefully. 

Due to the spectral gap 5 between and the rest of the spectrum of L p , we can fit 
a contour C around the origin in the resolvent set. Then (I4.17P shows that the spectrum 
cannot cross C for small k. A convenient choice for C is the rectangle 

C = (5 - e + i[-R, R]) U ([-R, 5 - e] + iR) U (-R + i[-R, R]) U ([-R, 5-e]-iR), 
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with R fixed independent of e, but sufficiently large. By Lemma El 



(4.18) sup (L ,p - z)' 



< oo. 



Expanding the resolvent of L kp in a Neumann series, 

oo 

(4.19) (L k , p - z)' 1 = ^(Z ,p - z)' 1 [(^o,p - Vp)(£o,p - z)~ 

n=0 

and using (I4.17P and (I4.18p . we see that there is r > such that if |k| < r, then C is in 
the resolvent set of L kp . However, the spectrum is a subset of the numerical range and the 
numerical range of L kp is contained in the set 

(4.20) {x + iy : x > k \y\ < C + 72} , 
with C = 2|| /i||oo + 2||V||. We conclude that 

(4.21) a(L k>p ) = S U Ei 

with E inside C and Si C {z : Rez > 5 — e}. 

It remains to show that So consists of a non- degenerate eigenvalue and to derive (14.141) 
and (I4.15p . For this purpose, consider the (non-Hermitian) Riesz projection 

(4.22) Q kiP = -L / -^^dz. 

The rank of Q^ p is constant so long as C remains in the resolvent set. Thus, Q^ p is rank 
one for |k| < r and So = {E p (k)} with associated normalized eigenvector $ kp in the one- 
dimensional range of Q k)P . Then, E p (0) = and $o,p = So®l. By the Feynman-Hellman 
formula, 

(4.23) diE p (k) = (^ Kp ,diL Kp ^ Kp ), 

from which it follows that VE p (0) = since VL k>p = iVfT k . p is off-diagonal in the position 
basis on % . Similarly, 

(4.24) 

didjEpQt) = ($k,p,^^Z k ,p$k, P ) + ($k,p,Qk^Z k ,p(^p(k) -Z ki p) _1 (l - Q k ,p)^L k)P $ kiP ) 

+ ($ k ,p,Q k ^L kiP (^p(k) - L Kp y l (l - Q k , p )9 i L kiP $ k) p) 

The first term on the r.h.s. vanishes at k = and the remaining two terms give (j4.15p . 
Because the form on the r.h.s of (14.151) is positive definite, the non- degeneracy condition on 
the hopping (Assumption 3) gives that didjE p (0) is positive definite. □ 

It follows from Lemma [3] and the sectoriality (12.31) of B that the semigroup e~* Lk,p satisfies 
exponential bounds (see [U Lemma 4]): 
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Lemma 4. Given e > there is C e < oo such that if k is sufficiently small, then 



(4.25) 



e -^,P(i-g kp ) 

for all p ; where Qk,p t/ie ranA; one Riesz projection (I4.22p onto £/je non- degenerate eigen- 
vector of Lk iP TOf/i eigenvalue near 0. 

4.2. Proof of Theorem [IJ, As in [4], it suffices to prove the theorem for p satisfying 
(4.26) ^2\po(x,y)\ < oo, 



since any initial density matrix can be approximated in trace norm arbitrarily well using 
such p . Assuming (I4.26p . note that 



(4.27) 



Po ; k,p(a;) = ^2 p (x- Nrj-y,-y)e 



ip-(:r— 7V77)— ik-j/ 



is uniformly bounded in £ 2 (A) as p varies through the torus: 



< 



.r,y 



(27r) d 

N d 
(2vr) rf 

N d 



dp(5 ® l,e- r "w, P p i) 



(4.28) Yl \Po-Xp( x ) 

By ( 13. 7ft . we have 
(4.29) 

^ e _1 ^ k ' x E (p Tt (x, x)) 

x 

(4.30) 

(431) ' 

By Lemma HI the integrand in ( I4.3ip is exponentially small in the large r limit 



f dp e -rtE p(WV -r ){So g, ^ g p ® 1) 



dp (6q <g> l.e-^/^pfl - Q i k;p )p 0; i k)P ® 1). 



(4.32) (So s 1, (1 - ^ p )e- rt ^ .u |p ® 1) 



< 



l - Q. k >- rt % 



P ;^k,p ® 1 



< ae- rt(5 - e) ->• 0. 



Regarding (14.301) . we have by Taylor's formula, 
(4.33) £ P (k/v^) = ^-J2 9 ^ E P^ k ^ 



since E p (0) = V£ p (0) = 0. Thus 
(4.34) e -T*£p( k /v^) 



o(l), 



12 



and 
(4.35) 

[ Z7T ) / nrd v T 

d p e 4E, J 3,¥p(o)k,k J ^ Oip(0) 

d 

JV 

since Qk p <5o £g> 1 — > 5o ® 1 as k — >■ and po;k, P (0) is continuous as a function of k. Letting 
Dij(p) = ~didjE p (0) and m(p) = j^dPo-,o, p (0) gives (12.91) and completes the proof. □ 
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